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Abstract 

Using computer simulations and scaling ideas, we study one-dimensional models of diffusion, 
aggregation and detachment of particles from islands in the post-deposition regime, i. e. without 
flux. The diffusion of isolated particles takes place with unit rate, aggregation occurs immediately 
upon contact with another particle or island, and detachment from an island occurs with rate e = 
exp(—E/kT), where E is the related energy barrier. In the partially reversible model, dissociation is 
limited to islands of size larger than a critical value i, while in the completely reversible model there 
is no restriction to that process (infinite i). Extending previous simulation results for the completely 
reversible case, we observe that a peaked island size distribution in the intermediate time regime, 
in which the mean island size is increasing, crosses over to the theoretically predicted exponentially 
decreasing distribution at long times. It contrasts with the partially reversible model, in which 
peaked distributions are obtained until the long time frozen state, which is attained with a crossover 

.3 

time r ~ — . The mean island size at saturation varies as S sa t ~ 2i + Ce (C constant), while the 
completely reversible case shows an Ahrrenius dependence of the mean island size, S ~ e -1 / 2 . 
Thus, for different coverages, the effect of the critical size i on the geometric features is much 
stronger than that of e, which may be used to infer the relevance of size-dependent detachment 
rates in real systems and other models. 

PACS numbers: 68.43.Jk, 81.10.Aj, 68.55.Ac 
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I. INTRODUCTION 



The early stages of formation of a thin film (submonolayer regime) is dominated by the 
deposition and diffusion of adatoms and the formation and growth of islands, through nucle- 
ation, coalescence and capture of new adatoms. This problem has long been of fundamental 
and technological interest [l| because those processes also have a strong influence on the 
later stages of growth. More recently, this interest increased for the possibility of formation 
of novel nano-structures. 

Systems where diffusion-mediated nucleation of islands competes with deposition of par- 
ticles have been extensively studied both experimentally and theoretically [J Q] . From 
the theoretical point of view, island formation (in one and two dimensions) has mainly been 
treated using scaling theories, rate equations and simulations. The models of deposition, 
diffusion and aggregation may be separated in three groups according to the conditions for 
the islands stability: models with completely irreversible aggregation 5), (], Q], with com- 
pletely reversible aggregation 8|, [9|], or those which assume the existence of a critical island 



The critical island size 



size i, above which islands are stable against dissociation 
is usually expected to represent thermal effects and the substrate geometry; for instance, 
i = 2 corresponds to stable trimers in a triangular lattice and % — 3 to stable tetramers in 
a square lattice |icl |. Other mechanisms that make island growth to be size-dependent may 
also justify such conditions. 

Another important problem is the nucleation and growth of islands in the post-deposition 
regime, i. e. after the particle flux has stopped. Diffusion and nucleation in this regime are 



observed in several experiments 



13 



14J ] . One possible situation is growth at low enough 



temperatures so that adatom diffusion is negligible during deposition, but with adatom 
rearrangement into islands occurring during a much longer time interval after the flux has 
been shut off Q]- This is certainly relevant for technological applications in which a 
system needs to have its properties preserved for long times. Other possible situation is that 
in which the sample is annealed pj5j. For the case of completely irreversible aggregation 
, the post-deposition regime was studied in one and two dimensions by several authors 



1) 



16,[l7,[l^]. The case i = 3 was also considered in two- dimensions by Tataru et al [it}]. Only 
for a point island model [3] , the case of partially reversible aggregation in a one-dimensional 
substrate has been treated in the post-deposition regime for several values of the critical size 
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1, 2 and 3). On the other hand, but also in one dimension, this post-deposition regime 



has been studied for completely reversible aggregation with an extended island model 20] . 

In this work, we consider one- dimensional models with critical size i for islands coarsening 
in the post-deposition regime, focusing on the differences from the completely reversible 
situation (i = oo), whose results from previous works are extended. First, the particles are 
randomly deposited on a one- dimensional lattice, and diffusion, aggregation and detachment 
mechanisms subsequently take place. Diffusion of isolated particles occurs with rate D = 1 
(probability 1/2 to move to each neighboring empty site), aggregation of a free particle 
immediately takes place upon contact with another particle or cluster, and detachment 
occurs with rate e if the cluster size does not exceed the critical size i (Fig. lb). The 
unit rate for hopping of isolated particles defines the time scale of the problem, and the 
detachment rate e is of order exp(—E /ksT), where E is an energy barrier for dissociation. 
The time evolution of island size will be analyzed for different detachment rates, coverages 
and critical sizes, and the island size distributions in different regimes will be obtained. One 
of our aims is to know how suitable is the approximation of finite critical size to describe 
the island coarsening without particle flux in different time regimes. 

At first sight, the restricted dimensionality considered here seem to have limited interest, 
but there are experimentally important cases of aggregate confinement, such as nucleation 



processes occurring along surface steps 
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22j. Moreover, some experiments show atoms 



diffusing along steps without detaching from them at room temperature 23| , the formation of 



one-dimensional nano-structures on patterned substrates [24 . l2q and the formation of single- 
atom- wide metal rows due to inhibition of aggregation at their sides 26]. Consequently, 
exploring nontrivial differences between these relatively simple models may be important 
for comparison with experimental data, as well as to design more complex models. This is 
illustrated, for instance, in recent work on two-dimensional irreversible island growth j^ . 
For these reasons, we will perform a careful quantitative analysis of the island properties 
and the time scales in the models presented here. 

The rest of this work is organized as follows. In Sec. 2 we present a brief review of previous 
results for the fully reversible model and simulation data for the island size distribution, 
whose shape shows significant changes in time. In Sec. 3 we discuss our results for the 
model with partially reversible aggregation. In Sec. 4 we summarize our results and present 
our conclusions. 
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FIG. 1: Diffusion and detachment processes of (a) the completely reversible model and (b) the 
model with a critical size. Detachment from islands with sizes greater than % is forbidden (denoted 
by x). 

II. THE COMPLETELY REVERSIBLE MODEL 



For the completely reversible model in the post-deposition regime 20[], the particles are 
randomly deposited on a one-dimensional substrate at t = and immediately begin to 
diffuse according to the conditions of Fig. la, but with no restriction to detach from an 
island (i = oo). 

The coarsening in this system was shown to be separated in three main regimes [^(J: (i) 
early fast attachment of isolated particles to each other; (ii) an intermediate regime in which 
detachment sets in, the mean island size increasing as 



S ~ (et) 



1/3 



(1) 



in the limit of small e; (iii) a slow (diffusive) evolution of the average cluster size as 



S(t) = Soo - C/t 



1/2 



where 



S n 



-1/2 



(2) 



(3) 



for small e. This last regime is diffusion-limited in the sense that the distances between 
clusters are still large so that a detached particle probably returns and reattaches many 
times before effectively diffuse to the nearest cluster. Actually this system does not attain 
neither an equilibrium nor a steady state. The behavior in regime (iii) and the corresponding 
island size distribution were predicted by a solution of the master equation in an independent 
interval approximation, while the behevior in regime (ii) was derived using scaling arguments 

y. 
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Extending previous kinetic Monte Carlo simulation, we confirmed the theoretically pre- 
dicted island size distribution in the asymptotic state [regime (hi)] [20! ]: 



P(x) 



X 



(4) 



r ''" 2(1-6) * s ^ e avera S e s i ze f° r coverage 6. This 



where x is the island size and S ~ f -i/2_9_ 1 / 2 j 
exponentially decreasing distribution is shown in Fig. 2 for e = 1/128 and 9 = 0.5, at a 
sufficiently long time t = 2 x 10 5 where the average cluster size is S = 9.53 (S^ = 9.99 in 
this case). However, the distribution at early times [regime (ii)] shows a peak at a certain 
typical size S ~ 4, as illustrated in Fig. 2 for the same model parameters and t = 3000. 
This new result was not predicted by the analytical solution which explained regime (iii). 
On the other hand, it is probably the reason for the successful derivation of Eq. ([1]) with a 
scaling approach that assumed the existence of a typical cluster size in the system in regime 
(ii). This assumption fails in the case of a very large number of small clusters, such as the 
asymptotic limit (iii). 
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FIG. 2: Island size distributions for the completely reversible model, for e = 1/128, 8 = 0.5 and 
L = 8000. At time t = 3000, the distribution presents a peak (squares), and at time t = 2 x 10 5 
the distribution decreases monotonically (circles). Averages over 1000 realizations. 



The crossover from regime (ii) to regime (iii) is expected to take place at times of the 
order e -5 / 2 , which can be obtained by matching the typical sizes in Eqs. (CQ) and ([3]). The 
above results show that this is the expected time for a crossover from a peaked island 
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size distribution to a monotonically decreasing one. This difference may be important for 
comparisons with experimental data and, as will be shown below, for a test of the reliability 
of the assumption of a critical island size. 

At this point, the comparison with results for reversible models where deposition competes 
with aggregation is also important. Without accounting for the post-deposition behavior, 
peaked distributions are observed in several conditions [8] . Moreover, they were also obtained 



by Lam and co-workers 



16j in the post-deposition regime, although a systematic search for 



an asymptotic regime was not performed there. Since this type of distribution is obtained 
in regime (ii) of our model, one naturally raises the question whether such crossover is also 
present in two dimensions. On the other hand, monotonically decreasing size distributions 



261 ]. and were explained by a highly 



of one-dimensional islands were obtained in Ref. 
anisotropic model of diffusion and completely irreversible aggregation (i — 1) in a square 
lattice. 



In models with conserved density, it is interesting to mention the case of Ref. [28|j, m 
which increasing the concentration of amphiphiles leads to a crossover from monotonically 
decreasing cluster size distributions to peaked ones. Moreover, depending on the particular 
spin dynamics, both types of distributions were found with same temperature and concen- 
tration. 



III. THE PARTIALLY REVERSIBLE MODEL 

In this model, after random deposition of a coverage 9 at t = 0, diffusion, detachment and 
re-attachment of particles continue until all islands are stable (size larger than i), following 
the rules of Fig. lb. We simulated this model in lattices of length L = 8000, which are large 
enough to avoid finite-size effects, with coverages 9 = 0.1, 0.5 and 0.9. The critical sizes 
analyzed here range from 2 to 70 and the diffusion rates e range from 10 _1 to 10~ 4 . 

We observed a sequence of four regimes of aggregation in the simulations: (a) fast at- 
tachment of isolated particles to each other to form islands (b) an intermediate regime in 
which detachment sets in, allowing further coarsening, but the size of the islands are typi- 
cally below i, thus with negligible influence of the critical size; (c) a crossover region, which 
begins when a large quantity of stable islands are formed; (d) a frozen state where there is 
no more isolated atoms nor unstable islands. Regimes (a) and (b) are equivalent to regimes 



6 



(i) and (ii) of the completely reversible model, which is justified by simulation results. 

In all cases, the mean size S increases in time and attains a saturated value S sa t that 
increases with i, for fixed e and 9. This is illustrated in Fig. 3a, where we compare the time 
evolution of S for i — 10, i — 15 and % = 50, with e = 10~ 4 . In Fig. 3b we compare the 
result for i = 50 with that for the completely reversible model with the same e. 

The results of Fig. 3a show that small values of i have drastic effects at intermediate times. 
Due to the stability of small clusters, the mean size S increases much earlier and exceeds 
the mean size for large values of i. For instance, S for % = 10 is larger than S for i = 50 in 
a time interval which extends for more than one decade in log (t) . The very slow growth in 
regime (b) [or (ii) in the completely reversible case] contributes to this effect. Consequently, 
the short-time behavior of S does not reveal the true asymptotics (S sat increasing with i) in 
cases where the critical size is a realistic assumption. 

These results are related to the fact that the intermediate regime (b) is clear only for 
relatively large i, otherwise the crossover to saturation takes place immediately after for- 
mation of small islands. Moreover, an increase of the mean cluster size similar to the fully 
reversible model [Eq. (TJ regime (ii)] is observed only for very small e, typically e < 10~ 4 . 
This illustrated in Fig. 3b, where we show a line with slope 1/3, which is nearly parallel to 
the data for both models in that regime. 

Contrary to the long time diffusive behavior of the completely reversible model (Eq. |2J), in 
the present case the crossover to the frozen state, regime (c), shows the expected exponential 
decay to a saturation value, as 

S(t) = S sat -e^\ (5) 

This is illustrated in Fig. 4 for i = 5, e = 10~ 2 and 9 = 0.5. Smaller values of i show much 
longer and well defined convergences to the saturation sizes. 

For fixed i, the mean island size in the frozen state, S sat , increases approximately as 

S S at = S + Ce, (6) 

typically with small C < 100, which corresponds to a very weak dependence on e in the low 
temperature regime. This is illustrated in Fig. 5, where we show S versus e for i = 5 and 
different values of the coverage 9. Except for coverages very near 1, the saturation value for 
low temperatures (e <C 1) is nearly So ~ 1i\ for instance, So ~ 10 in Fig. 5 for 9 = 0.1 and 
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FIG. 3: Time evolution of the mean island size S, for = 0.5, e = 0.0001 and L = 8000 (a) 
for different values of the critical size: i = 10 (dashed) , i = 15 (dotted), averages over 1000 
realizations, and i = 50 (solid), averages over 50 realizations; (b)for i = 50 (solid) and for the fully 
reversible aggregation model (dashed). The dotted line has slope 1/3. 

9 = 0.5. This expected proportionality between So and i is also observed for larger values 
of i. 

This result contrasts to the remarkable dependence of the mean island size on e for the 
completely reversible case [regime (iii), Eq. [3]. Here, the critical island size is the main 
parameter to determine the features of the system in the saturation regime, while e has 
an important role in determining the time scale of the problem. It is also the parameter i 
which is the more closely connected to other geometric features, such as size distributions 
(discussed below). 

The crossover time r in Eq. ([5]) can be estimated from the previously discussed features 
of both models. First recall that Eq. ([I]) governs the time evolution of the mean island 
size in intermediate times when both models have similar features (small e, large i). On the 
other hand, at the crossover time r, this size must match the saturation value S ~ 2i for 
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FIG. 4: S sat - S(t) for % = 5, 9 = 0.5, e = 1(T 2 , L = 8000, averaged over 1000 realizations. 
S sat -S(t) = 7.12e^/^ 



The result r ~ 1/e is expected because the detachment rate defines the time scale of the 
model; it is confirmed by simulation for small e with good accuracy. However, much more 
interesting is the nontrivial dependence of r on j. In Fig. 6, it is confirmed by simulation 
with large values of i (i = 30 to i — 70), for e = 10~ 4 and 9 = 0.5. Notice that, despite the 
small range of values of i, i 3 varies by a factor larger than 5 in Fig. 6. This scaling certainly 
requires relatively large values of i and small e to be observed. On the other hand, this 
analysis reinforces the conclusion that the intermediate time regime of this model, regime 
(b), has negligible differences from the model with completely reversible aggregation. 

The final distribution of island sizes (frozen state) has a pronounced peak for small i, 
due to the stability of small clusters that stop feeding the largest ones at early times. As 
expected, the distribution becomes larger when i increases, as illustrated in Fig. 7. For fixed 
i and 8, the distribution becomes narrower for smaller e. 

In order to illustrate the remarkable difference from the fully reversible model at long 



the present model (S ~ i is certainly intuitive). Thus we obtain 



r ~ — 

e 




(7) 
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FIG. 5: Mean island size at the final state, S sa t, as a function of e, for i = 5, 9 = 0.1 (circles), 
9 = 0.5 (squares), 6> = 0.9 (triangles). 

times, we compared the island size distributions for the same e and S sat ~ S^, using 9 — 0.5. 
For instance, for i = 15 and very small e we obtain S sat ~ 32, and the same value for Soo is 
obtained in the fully reversible model with e = 0.00049. The distributions of both models 
are compared in Fig. 8a, with the same value of e and i — 15 for the model with critical 
size. The same comparison is also performed in Fig. 8b for i — 5, 9 — 0.1, and e = 0.000845, 
where S sa t ~ 9.5 in both models. Thus, from the shape of the size distribution at very long 
times, we are able to test the reliability of each model for a given system. 

IV. CONCLUSION 

We studied island coarsening in the post-deposition regime in one-dimensional models 
with completely and with partially reversible aggregation, the latter with critical size i. 
Extending previous work on the fully reversible case, here we showed that the island size 
distributions are very different in the regime of island coarsening and in the long time limit, 
crossing over from a peaked distribution to an exponentially decreasing one. For the model 
with critical size i, we systematically studied the influence of this parameter, the detachment 
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FIG. 6: Dependence of the crossover time r on the critical size i. For i = 30, 40, 50, 60 and 70, 
6 = 0.5, L = 8000, averaged over 50 realizations. 

rate e and the coverage 9 on the mean island size and the island size distribution . Contrary 
to the long time diffusive increase of the average island size in the fully reversible model, 
for the partially reversible model the crossover to a frozen state shows an exponential decay 

.3 

to the saturation value, with crossover time r ~ — , a relation which follows from the fact 
that, at intermediate times, both models show the same coarsening. For fixed i, the mean 
island size in the frozen state has a very weak dependence on e, in contrast to the completely 
reversible model, where it varies as S ~ e -1 / 2 . Another particularly interesting result is that 
the peaked island size distribution at intermediate times for both models crosses over to a 
long time exponentially decreasing distribution only for the completely reversible model. 

We believe that some of these features may motivate comparisons with real systems 
with effective one-dimensional behavior (see e. g. the recent discussion in Ref. 29j. For 
instance, the completely reversible model predicts an Ahrrenius temperature dependence of 
the mean island size (Eq. [3] with e ~ exp (—E/ksT)) in the long time limit, i. e., when that 
size does not show significant time increase. Deviation from this dependence may suggest 
some size-dependence of detachment rates, and would eventually justify the assumption 
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FIG. 7: Island size distribution for 9 = 0.5, e = 0.1 i = l(circles), 2 (squares), 5 (diamonds), 
10 (triangles), 15 (stars). 



of a critical island size or another restriction to the growth of large islands for modeling 
the island coarsening. Moreover, peaked island size distributions may indicate deviations 
from the long time conditions of the completely reversible model, which may be interpreted 
as a finite-time behavior [regime (ii) of the same model] or as an effect of size-dependent 
detachment rates. 

There are also many recent examples in which the solution of simple statistical models in 
one-dimension anticipate features of the more realistic two- or three-dimensional systems, 
such as non-equilibrium multilayer growth models [30|], magnetic systems with rough surfaces 



3l| and other aggregation models 32] . Thus, further systematic studies of island coarsening 



in the post-deposition regime in two dimensions are certainly motivated by the present work. 
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FIG. 8: Island size distributions for the same values of S sa t and e, using (a) 8 = 0.5 and e = 0.0049, 
S sat ~ 32 , for i = 15 (circles) at the final state, and for the fully reversible model (triangles), at 
t = 5 x 10 7 (asymptotic regime), (b) 9 = 0.1 and e = 0.000845, S^t as 9.5, for i = 5 (circles) at 
the final state and for the totally reversible model (triangles), at t = 3 x 10 7 (asymptotic regime). 
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